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ABSTRACT 

We have searched for nonastrophysical emission lines in the optical spectra of 
577 nearby F, G, K, and M main-sequence stars. Emission lines of astrophysical 
origin would also have been detected, such as from a time-variable chromosphere 
or infalling comets. We examined ~20 spectra per star obtained during four 
years with the Keck/HIRES spectrometer at a resolution of 5 km s“^ , with a 
detection threshold 3% of the continuum flux level. We searched each spectrum 
from 4000 A - 5000 A for emission lines having widths too narrow to be natural 
from the host star, as well as for lines broadened by astrophysical mechanisms. 
We would have detected lasers that emit a power, P >50 kW, for a typical beam 
width of ~0.01 arcsec (diffraction-limit from a 10-m aperture) if directed toward 
Earth from the star. No lines consistent with laser emission were found. 

Subject headings: astrobiology, extraterrestrial intelligence — techniques: spec¬ 
troscopic 


1. Introduction 

Most searches for signals from extraterrestrial intelligent civilizations (SETI) have con¬ 
centrated on radio wavelengths (e.g. Tarter 2001, Gullers 2000, Werthimer et ah 2000, Leigh 
and Horowitz 2000). However, a preferred wavelength for SETI work has proved elusive. The 
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motivation for this project stemmed from the suggestion of Schwartz and Townes (1961) that 
advanced civilizations might use optical or near-infrared lasers as an efficient means of inter¬ 
stellar communication (Townes 1998, Lampton 2000). The advantages of optical wavelengths 
over radio include the higher bandwidth and the narrower diffraction-limited beam size to 
promote private communication. In addition, the technology needed to produce lasers with 
detectable powers is not far from current Earth technology (§2.3). 

Optical SETI (OSETI) projects are currently designed to detect highly energetic, pulsed 
lasers and continuous-wave lasers. Such efforts include dedicated searches, piggyback searches, 
and data mining. Many OSETI projects are designed to detect pulsed lasers using fast pho¬ 
tometry (Howard et ah 2000, Bhathal 2000, Werthimer et ah 2001). Two or three photode¬ 
tectors are employed to search for laser pulses having a duration of ~1 ns in broadband 
visible light. A high-powered nanosecond pulsed laser will outshine its host star during its 
brief pulse because at most one photon will be received (for a typical 1-m class telescope) 
from the star during the 1 ns duration. 

Alternately, interstellar optical communication may occur via nearly continuous laser 
emission. In that case, lasers of even moderate power, tens of kilowatts, can outshine the 
host star when focussed by large apertures during transmission (see §2.3). Here we report a 
search for ultra narrow emission lines such as predicted by the Schwartz-Townes suggestion. 


2. Optical SETI with High Resolution Spectroscopy 

We examine high resolution spectra of 577 F, G, K, and M stars, most of which reside 
within 50 pc. We have taken ~12,500 spectra using the HIRES echelle spectrometer on the 
10-m Keck 1 telescope. All target stars are listed in Nidever et ah (2002). The stellar spectra 
were obtained during the period 1997-2001 in a search for extrasolar planets using relative 
radial velocities. We are indebted to Drs. R. P. Butler and S. S. Vogt for their contributions 
in obtaining these Keck spectra. The multiple spectra from 4000 A - 5000 A are compared 
to detect sharp differences such as those caused by nonnatural emission lines, as described 
in the following subsections. 


2.1. High Resolution Echelle Spectra 

The acquired echelle spectra contain 34 spectral orders, covering wavelengths 3800 A- 
6100 A with resolution A/AA = 70,000 (see Vogt et ah 1994, Vogt et ah 2000, Butler et 
ah 2000). Not all orders are used here because the starlight is sent through a glass cell 
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filled with iodine vapor before entering the spectrometer for use in setting the wavelength 
scale in the planet search. The iodine lines contaminate the stellar spectrum at wavelengths 
greater than 4980 A. The remaining usable orders, 72 through 89, contain 1000 A ranging 
from 3980 A - 4980 A. This region includes Balmer /3 and 7 , but not the Call H&K lines. 
The pixel width is given by A/AA = 145,000 , corresponding to 2.1 km s“^ per pixel or 
~0.032 A per pixel. The actual spectral resolution is only 70,000 because the instrumental 
prohle has a width of 2.3 pixels. 

The recorded two-dimensional raw spectra are reduced in the usual way to one-dimensional 
spectra. The bias is subtracted from all images including both the stellar spectra and the 
flat-held exposures. Each spectral order recorded by the CCD has a width proportional to 
the diameter of the stellar seeing disk along the length of the slit. The typical seeing at 
Keck/HIRES of 0.7-1.2 arcsec yields orders of width ~4 CCD rows (each “row” being the 
on-chip sum of neighboring original CCD rows). Cosmic rays are removed by searching for 
7-a sharp intensity spikes at one row compared to the intensity in a neighboring row (at the 
same wavelength). We are indebted to Jason Wright for this cosmic-ray removal algorithm. 
The stellar spectra are divided by the flat-held exposures. The two-dimensional spectral 
orders are compressed to one dimension by summing the recorded counts within the ~10 
rows centered on each order. We thank Eric Nielsen for the careful reduction of CCD images. 
Figure 1 shows the results of this process. The signal-to-noise ratio per pixel in orders 72 
through 89 of a typical 1-D spectrum is ~200. 


2.2. Overview of the Detection Scheme 

We search for laser emission lines in the reduced spectra by comparing each spectrum 
to a reference spectrum of the same star, searching for statistically signihcant diherences. 
Over time, a star’s spectrum should remain nearly constant except for minor huctuations 
due to photon noise and photospheric changes. A large diherence at one wavelength could 
signify the presence of a laser signal that had been detected during one observation but not 
during the other. 

We would fail to detect an emission line if it occured at the same wavelength in both the 
observed and reference spectra. A laser hred continuously from the rest frame of the stellar 
photosphere would go undetected. However, such a continuous laser would be detected in 
both spectra at different wavelengths if the laser’s velocity vector were to change between 
observations, e.g. due to orbital motion around its host star. The resulting Doppler shift 
permit detections in the difference spectrum. 
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For each star, a series of algorithms is applied to all of its spectra. First a “reference” 
spectrum is chosen. Then, one by one, the other “test” spectra are aligned in wavelength with 
the reference for comparison (see §3.2). Candidate laser lines are found in the test spectra 
by subtracting the reference spectrum and searching for positive differences, or “spikes.” 
One test spectrum is then used as a secondary reference to check for spikes in the primary 
reference spectrum. The spikes are subjected to a selection criterion based on their width, 
to ensure that their shapes are consistent with the instrumental prohle of HIRES (see §3.3). 
Spikes that are narrower than the instrumental prohle are likely to be caused by cosmic rays 
hitting the CCD that eluded detection in the CCD reduction process described in §2.1. 

A monochromatic laser line would have an intrinsic width narrower than the resolution 
of the spectrometer. However, as the laser light passes through the HIRES optical system, 
it would be smeared out by the HIRES point spread function (PSF), giving the detected 
emission line an apparent distribution in wavelength characteristic of the HIRES PSF. Thus 
candidate laser lines can be distinguished from most cosmic rays, most CCD haws, and 
stellar emission lines by their predictable spread in wavelength that would be merely the 
HIRES PSF, as described in detail by Valenti et al. (1995). Notably, emission lines having 
widths that are narrower than possible from a star (with its thermal, collisional, turbulent, 
and rotational broadening mechanisms) would stand out as possibly “nonastrophysical” in 
origin. All of the resulting candidate emission lines of any width, excluding the apparent 
cosmic rays hits, from the spectra of each star are stored for further analysis as described in 
§4. In sections §3.1, §3.2, and §3.3, we describe in detail the processes by which we searched 
for emission lines. 


2.3. Detection Limits for Laser Emission 

To be detected, a laser signal must compete against the light of its host star. A laser 
beam carries three advantages over starlight. First, the laser light is concentrated in a 
narrow cone of tiny solid angle relative to starlight, which is omnidirectional. Second, a 
monochromatic laser line falls into one resolution element, with a width of ~2.3 pixels, while 
starlight, composed of a wide range (~4000 A) of wavelengths, is spread out into roughly 
50,000 resolution elements. Lastly, the laser needs only to outshine the photon noise of 
the stellar flux, not the flux itself, within that resolution element since we are looking at 
differences in the spectra. To estimate these advantages of the laser signal, we consider the 
solid angles of the starlight and a diffraction-limited laser beam. 

f^star 


47r 


( 1 ) 



-5- 




laser 



( 2 ) 


Here A is the wavelength of the laser light and a is the aperture size of the laser transmitter. 
The fraction of starlight contained within the laser beam is then given by 

We note that the intrinsic width of a laser line will be much narrower than the HIRES 
resolution of ~0.075 A, thereby placing most all of the laser light into a single resolution 
element. The stellar spectral energy distribution may be approximated by a blackbody curve. 
For a Solar-type star, the SED has a width of ~4000 A, yielding the number of resolution 
elements, Wtar, into which the stellar luminosity falls: 


4000 A 
0.075 A 


5.3 X 10^ 


(4) 


The fraction of starlight contained within one resolution element is therefore given by 


A 
a 

The laser light needs only to outshine the stellar photon noise since we are detecting 
differences in a star’s spectrum over time. The photon noise relative to the flux itself, a, 
for one spectrum is given by 1/(S/N) where S/N is the signal to noise ratio per resolution 
element. Since the laser must compete with the photon noise of two spectra in one resolution 
element, we have 
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Finally we hnd that the laser power necessary for an na detection is 


1 /Ay 1 nV2 

~^\a) 5.3 X 104 S/N 


( 6 ) 


(7) 


where Tstar is the luminosity of the star. As a representative case, we take A=4500 A, 
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S/N=200, n=6, Lgtar = -^^sun = 3.8 X 10^® W, and assume Keck to Keck communication so 
a=10 m, which gives 


Plaser = 50 kW (8) 

Thus for this representative case, one can detect a 50 kW laser. 

The previous calculation assumes that the laser is on continuously. If, however, a 
nanosecond pulsed laser signal were transmitted. Equations 7 and 8 would need to be mul¬ 
tiplied by the factor jv ,^xio -9 s ''^^ere Atexp is the exposure time of the observation and 
TVpuises is the number of pulses received in that time. Using a typical exposure time of 10 
minutes and assuming 1 pulse is received, the minimum necessary power of a detectable 
nanosecond pulsed laser is Ppuised = 3 x 10^® W. Thus, the minimum necessary energy of a 
detectable nanosecond pulsed laser is 30 MJ. 

For comparison, the most powerful continuous laser on Earth is the Free Electron Laser 
(EEL) developed at Jefferson Laboratory in Virginia. The EEL has reached an average 
power of 1.7 kW and the laser is expected to achieve average powers up to ~10 kW with 
an upgrade (Jefferson Lab 2001). Lawrence Livermore National Laboratory produced the 
world’s most powerful pulsed laser, the Petawatt laser. During its 440 femtosecond pulse, it 
delivered a power greater than 10^® W (LLNL 2001). 

We can also determine a minimum detectable flux, Tiaser, 

^ (9) 

-^laser 

where Aiaser is the area of the laser beam at the Earth and given by 

Aaser = ^Q) (10) 

where d is the distance from the laser to Earth. The flux limit can now be expressed as 


4 1 71^/2 Lstar 

7r5.3 X 10^ S/N ATid'^’’ 


( 11 ) 


where Lst^.T./Airdi^ is the flux of the star. If we take d = 100 light years with the previously 
adopted values, the minimum flux of a detectable laser is 
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Fiaser = 3.4 X lO"'" W/m^ (12) 

Thus, we can detect monochromatic laser light from nearby stars if the flux is at least 
3.4 X 10“^’^ W/m^ at Earth. 


3. Spectral Analysis 

A detailed description of the algorithms developed to analyze the spectra is given in 
this section. 


3.1. Choosing Reference Spectra 

Two reference spectra are required for each star. The primary reference spectrum is 
compared to all the other test spectra. The second is compared to the primary reference 
spectrum. Since only positive differences are considered, any candidate laser lines in the 
hrst reference would not be flagged during its comparison with the test spectra. A second 
reference, one of the previous test spectra, is used to account for this possibility. 

The two reference spectra are chosen as follows. Almost every star has at least one 
“template” spectrum. In the planet-search project, the stellar radial velocities are measured 
with respect to the template spectra. These observations are of good quality with high 
signal-to-noise ratios. If there are multiple-template spectra for a star, the two with the 
highest signal to noise are chosen as the reference spectra. If there is only one template, it 
is used as the primary reference spectrum and another observation with high signal to noise 
is chosen as the second reference. Some stars which are new to the extrasolar planet search 
catalog do not yet have templates. For these stars, the two observations with the highest 
signal to noise are chosen as the references. 


3.2. Aligning Spectra 

Aligning a pair of spectra precisely in both wavelength and flux is important. If the 
spectra are exactly the same, except for minor fluctuations due to photon noise, the fractional 
standard deviation, a, of the difference spectrum is expected to be 




where (S/N)i and (S/N )2 are the signal-to-noise ratios of the two spectra. For two spectra 
with signal-to-noise ratios of 300, a can, in principle, be as low as 0.005. It is important to 
align the two spectra well enongh snch that the observed error is as close to the theoretical 
error as possible in order to detect the smallest differences. 

We line np a pair of spectra order by order. Various normalization and shifting algo¬ 
rithms are applied to both the reference and test spectra, as follows. The hrst step is to put 
the orders on roughly the same flux scale by dividing by the continuum. The continuum is 
found by first partitioning an order into 10 bins. Each bin is assigned a temporary continuum 
at the 85th percentile in flux. A third-order polynomial is £t to these individual continua 
to normalize the entire spectral order. 

Each spectrum is Doppler-shifted to correct for the Earth’s orbital motion, which causes 
a maximum wavelength shift of ~1 A corresponding to ~30 pixels. We determine the 
Doppler shift of the test spectrum relative to the reference spectrum by minimizing and 
then shift the test spectrum so that the two are aligned. The final shifts are precise to 
~0.01 pixel, and fractional pixel displacements of the spectra are accomplished with a spline 
interpolation. This shifting is first carried out on segments of length 280 pixels, and then 
it is done on a 40-pixel scale (“chunks”). If any spectrum contains cosmic rays (or laser 
emission lines), it is difficult to align. Such features are removed temporarily (replaced with 
a spline interpolation) to carry out the alignment. Operating on chunks of 40 pixels allows 
us to account for the wavelength dependence of the Doppler shift and to correct any slopes 
in the flux with wavelength. 

The resulting (temporarily spike-free) reference and test spectra are subtracted to yield 
a difference spectrum, from which the standard deviation is determined and recorded as the 
“noise” for each spectral chunk. Finally the reference and test spectra are subtracted with all 
spikes included, leaving a difference chunk which can be searched for laser-line candidates. 
Figure 2 shows a chunk of difference spectrum both with (top) and without (bottom) flux 
spikes of unidentified origin. 

3.3. Criteria for Laser Line Candidates 

A laser beam would travel through the telescope and the HIRES spectrometer optics just 
as starlight would. Moreover, the laser beam travels through the same column of Earth’s 
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atmosphere, acquiring the same wavefront distortions and seeing prohle as does the star. 
Therefore, an emission line from a laser would span the whole width of a spectral order 
(perpendicular to dispersion) as would the stellar spectrum on the CCD. While the intrinsic 
width of the laser line in wavelength would be smaller than the resolntion of HIRES, it would 
be broadened by the point-spread fnnction of HIRES. The B1 slit used in HIRES projects 
to a width of 2.0 pixels on the CCD. The instrnmental prohle of HIRES from both its optics 
and the slit width is known from carefnl measnrements (Valenti et al. 1995) and from simple 
Gaussian htting of emission lines from a thorinm-argon lamp. The instrnmental prohle with 
our HIRES setup has a width of 2.3 pixels (FWHM). Thus, we expect any laser line to have 
an observed width of 2.3 pixels. 

We search for candidate emission lines in the stellar spectra as follows. Each 40-pixel 
chnnk of diherence spectrnm has a corresponding fractional photon noise, a, associated with 
it as described in §3.2. We pass throngh all the chunks and note spikes having heights above 
6(7. This threshold is chosen to exclnde essentially all noise bnt detect a signal from a laser 
with the lowest detectable power. 

We also reqnire that candidate laser lines exhibit a width that is at least as broad 
as the HIRES instrnmental prohle, with EWHM=2.3 pixels, allowing us to discriminate 
against (i.e., reject) cosmic-ray hits. Most cosmic rays prodnce electrons in a single pixel or 
occasionally in two neighboring pixels, clearly narrower than the HIRES PSE. To implement 
this “width” criterion, we reqnire that the shorter of the two adjacent pixels on either side 
of each G-a spike contain at least 30% of the photons in that central spike. Examination of 
the known HIRES PSE and thorinm comparison lines (see Fignre 3) shows that this width 
criterion wonld be met by almost any emission line that had been broadened by the HIRES 
PSE. Any spike in the difference spectrnm that does not meet both the G-a height criterion 
and this width criterion is rejected from fntnre consideration as they are probably cosmic 
rays. Spikes that do meet both criteria are stored in a database that inclndes the location 
and characteristics of that spike, for later statistical analysis (§4). We note that emission 
lines having widths greater than the PSE will be detected and stored as candidate signals. 
Thus our survey would detect intermittent emission lines of any origin, laser or astrophysical, 
if they are sufficiently strong. 


4. Statistical Analysis of Candidates 

The spectral analysis described in §3 yields a database of surviving candidate emission 
lines in each of the ~12,500 Keck spectra of onr 577 target stars. A typical star has ~13 
surviving candidate emission lines in its ~20 spectra. For each candidate, the database 
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contains the star, date of observation, the wavelength, the peak number of photons, and 
the width of the candidate line. Detailed examination shows that most (if not all) of the 
surviving candidate emission lines are simply cosmic ray hits that we had failed to reject. 

To discriminate a true emission line from these cosmic rays, we apply a statistical 
analysis on the collection of emission-line candidates from all spectra (typically 20) for a 
given star. We search for an overabundance of candidates (“hits”) in the vicinity of any 
particular wavelength. In contrast, cosmic rays would hit at random locations on the CCD, 
and would not favor any particular wavelength. A clustering of hits near one wavelength 
from all the spectra of one star could indicate laser emission. 

We note that the Earth’s orbital motion will not cause a Doppler shift of the laser line 
in our analysis, as we have already forceably shifted the reference and test stellar spectra 
onto the same wavelength scale (§3). However, acceleration of the emitter relative to the 
star (such as due to orbital motion) during the typical 4 years of observations would cause 
the laser line to be Doppler-shifted to neighboring wavelengths. 

We make a grid representing the CCD with the spectral orders partitioned into approx¬ 
imately 30 segments of roughly 100 pixels each. The segments overlap one another by half 
of their size to account for the possibility of an abundance of hits split into two neighboring 
segments. We record the number of emission-line candidates in each segment. 

For each star, we plot the number of 100-pixel segments against the number of hits in the 
segments. Most segments have no hits at all, and only a few have one hit. These infrequent 
hits (none or one hit out of ~20 observations in a wavelength segment) are consistent with 
random cosmic ray hits. To this histogram, we £t a Poisson distribution (see Figure 4), 
and determine the probabilities of finding an excess number of hits, n, in any wavelength 
segment. A star exhibiting an excess number of hits in a segment, i.e. a low probability of 
the events occurring randomly, would stimulate further investigation as a potential emission 
line. 


The fit of a Poisson distribution to the histogram of hits is done as follows. We fit the 
low-n (0, 1, 2) portion of the histogram and neglect the high-n tail. This approach avoids 
contaminating the Poisson mean, /i, with actual emission lines in the high-n tail. The first 
empty bin in the histogram, no, is determined and an estimate of the expected mean number 
of hits per segment, /i, is calculated neglecting those segments with n > Uq. 

The Poisson probability of finding n cosmic ray hits in any given wavelength segment 
is given by (Taylor 1982): 



11 


PM = (14) 

Using the estimated /i, we calculate Pfj,{n) for n = 0,1, ...,no — 1, and check if any of the 
probabilities are less than 1%. If so, we exclude those bins from the £t as well if no > 2. 
The fit is then applied to bins from n = 0 to some nmax hits. 

The estimated mean number of hits per segment may not necessarily yield a Poisson 
distribution that gives the best fit to the observed histogram. The best Poisson fit may be 
better determined by finding the lowest given by 
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where is the observed number of segments containing n candidate hits and E^ is the 
expected number with n hits, based on a Poisson distribution. The only free parameter is 
/i, the mean number of hits per segment. We vary fi to find the value corresponding to the 
minimum of The resulting Poisson probabilities are multiplied by the total number of 
segments and adopted as the best fit to the observed histogram. 

We now determine the probabilities of finding n hits in any given segment, extrapolating 
beyond nmax- Segments having an excessive number of hits with probabilities less than 0.01% 
are noted as serious candidate emission lines, as they occur in excess numbers above the 
background cosmic ray hits. Since there are about 500 wavelength segments on our grid of 
the CCD for each star, we expect to flag one in 20 stars spuriously, i.e. due to a chance 
coincidence of cosmic rays. 

We find that indeed 45 stars were flagged as having a wavelength segment with an excess 
number of hits at the 0.01% false-alarm level. Each of the surviving flagged candidate laser 
lines, of which there were 187, were examined by eye in the raw CCD images. 


5. An Apparent Emission—Line Detection 

Our search for intermittent emission lines revealed an apparent detection for the star 
HD29528 (KOV). Three separate spectra obtained on different dates contained strong can¬ 
didates at the same wavelength. The hits were noted in both of two overlapping spectral 
segments for this star. The probability that these three apparent emission lines were merely 
spurious cosmic rays occurring within a segment of length ~200 km s“^ is easily calculated 
from the Poisson occurrence of cosmic rays (§4). 



Figure 4 shows the observed histogram of the number of segments versus number of hits 
along with the best htting Poisson distribution. The probability of randomly obtaining 3 
hits in a segment for this star is 1.23 x 10“^ or 1 in 8 million. 

Figure 5 shows these candidate laser lines in the rednced spectra of the three observa¬ 
tions. The aligned reference and test spectra are shown and twice the difference spectra are 
below them. The candidate laser lines are in spectral order 83 and the central peaks are in 
pixel 1436 with respect to the reference spectrnm. 

We examined the raw spectra to see if the candidate laser lines span the whole width 
of the spectral order, indicating whether the light passed through the telescope and spec¬ 
trometer. Figure 6 shows the raw CCD image of the spectrum of HD29528 obtained on 
JD 2451793.1, centered on the hrst candidate laser line. The reference spectrnm is shown 
below it. Beneath the reference is the raw difference spectrnm making it easier to see the 
apparent laser line. The other two candidate laser lines fonnd in observations obtained on 
JD 2451882.98 and on JD 2451884.08 are not shown, bnt those apparent emission lines also 
span the width of the spectral order. Because all three hits occur at the same pixel (and 
hence wavelength) with respect to the star’s reference spectrnm, and because they span the 
entire width of the spectral order, the hits are unlikely to be cosmic rays. These are surviving 
laser line candidates thus far in the analysis. 

To further investigate these candidate emission lines in HD29528, we compared the 
spectrum of HD29528 to the solar spectrum using “The Solar Spectrum 2935 A to 8770 A” 
(Moore et al. 1966). The rest wavelengths of the lines are 4307.31 A and there are no listed 
absorption lines dne to the Earth’s atmosphere nor anomalous Solar emissions lines at this 
wavelength. 

We further examined the spectrnm of a rapidly rotating B-type star taken on the same 
night as the reference spectrnm for HD29528. The B star spectrum should have no prominent 
spectral features. A strange drop in flux of 20% is seen in the spectrum, exactly at pixel 
1436 in order 83, precisely where the hits were fonnd with respect to the reference spectrnm 
of HD29528 (see Figure 7). This clearly spurious feature in the B star strongly snggests that 
onr candidate laser lines are merely artifacts cansed by a flaw in the CCD. 

Indeed, a close examination of the raw CCD image of the reference spectrnm of HD29528 
shows the flaw in the CCD (see Figure 8). This flaw is certainly the sonrce of the false 
detection. 

As a side note to this false detection of an emission line, one may calculate the necessary 
power of a laser that conld have prodnced a “signal” similar to that detected. In Eqn 
7 we assnme a transmitter apertnre, a=10 meters. The other variables are known. The 
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wavelength of the 13a “signal” is 4307 A, and the signal-to-noise ratio of the spectra is 
~120. The luminosity of a KO star is 0.42 Lsun (Carroll and Ostlie 1996). Using these values 
in equations 7, we hud Piaser = 68 kW. This would be the required power of a laser to produce 
the spurious signal in Figure 6. This test supports our expectation that lasers having power 
of 50 kW would have been detected. 


6. Results 

For each of the 45 stars that exhibited an excess of candidate emission lines (above 
the cosmic ray background noise) we examined the original raw CCD images from HIRES. 
Approximately 200 raw spectra were examined to follow up on the 45 flagged stars. Stellar 
emission lines were detected in the spectra of two flare stars, HIP5643 (M4.5) and HIP92403 
(M3.5). Both stars exhibited varying emission of the Balmer line H/5. 

Interestingly, many candidate emission lines were located along the steep slopes of ab¬ 
sorption lines in the reduced spectra. However, the corresponding emission spikes were not 
apparent in the raw CCD image. Such emission spikes are apparently neither cosmic rays 
nor emission lines. Instead, it is likely that these spikes located in absorption lines are caused 
simply by tiny changes in the instrumental prohle of HIRES from one observation to the 
next. The changes in the instrumental prohle could result from different atmospheric see¬ 
ing of the star image at the slit (making the slit effectively a different size) or from small 
focus changes in HIRES itself. The steep sides of absorption lines have local intensities that 
are very sensitive to the instrumental prohle. It is also possible that the spikes located in 
absorption lines are residual numerical artifacts from the spectral alignment process since 
interpolation routines normally give the largest errors on steep slopes. Thus these candidate 
emission spikes located within absorption lines were rejected as spurious. 

All other candidate emission lines were visually inspected and found to be inconsistent 
with actual emission lines. Some were multiple, surviving cosmic rays that occurred, by 
chance, near the same wavelength. Some emission spikes were caused by bad pixels or haws 
in the CCD. No lines consistent with laser emission were found in any of our stars. 


7. Discussion 

We would have detected time-varying emission lines from any of the 577 main sequence 
stars above a threshold of a few percent of the continuum hux. Emission lines of constant 
intensity would also have been detected had they Doppler-shifted relative to the star during 
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the ~4 yr of observations. Lines having intrinsic widths from arbitrarily narrow to tens 
of km s“^ would have been detected. No narrow lines consistent with laser emission were 
detected. However, naturally occurring emission lines were found in the spectra of two stars. 

Emission lines in the spectra of Solar-type stars are rare between 4000 A - 5000 A. 
Such lines can arise occasionally in the chromospheres or coronae of the most magnetically 
active stars. In some (~20%) F, G, K-type main sequence stars, the cores of the Balmer 
absorption lines are “hlled in” by a few percent due to chromospheric emission (Herbig 1985, 
Robinson, Cram, & Giampapa 1990). We would have detected temporal variation in that 
chromospheric hlling, if it exceeded a few percent of the continuum level. No such variations 
were detected in any of our F, G, K-type main sequence stars. This nondetection indicates 
that such chromospheric variation in the cores of the strong Balmer lines, notably H/5 and 
Hy, did not exceed a few percent of the continuum level in these stars. Flare stars, most 
of spectral type MO or later, more commonly exhibit emission of H/?. Indeed, we detected 
time-varying emission of H/5 in two M-type flare stars. 

Time-varying emission lines from stars may also arise from infalling volatile material 
and comets, as is observed for f3 Pic (see for example Beust et al. 1998). Our search would 
have detected such astrophysical emission lines if they appeared or varied by more than a 
few percent of the continuum flux of the star. All such astrophysical origins for emissions 
lines would yield line widths greater than our resolution of 5 km s“^ due to the usual thermal 
and turbulent broadening, as well as possible rotational and collisional broadening. Hence 
such lines would be resolved. Nonetheless, no emission lines due to these effects were found 
in the 577 stars at our threshold of a few percent of continuum flux. 

The entrance slit of the HIRES spectrometer is rectangular with dimensions 0.57 x 3.5 
arcsec. At the typical distance of our stars of ~50 pc, this slit corresponds to 30-175 AU 
at the star. Thus, sources of emission lines residing within 30 AU of our target stars would 
have been included in our spectra. Of course, the emission beam must be pointed toward 
the Keck telescope for detection. Extinction from interstellar dust is negligible toward these 
target stars due to their proximity (d <100 pc). 

The motivation for this search stemmed from the suggestion of Schwartz and Townes 
(1961) that advanced civilizations might communicate from one colony to another within a 
planetary system, or from one planetary system to another, by using optical lasers as the 
carrier. The advantages of optical communication over radio include its higher bandwidth 
and its smaller diffraction beam to promote private communication. 

A nondetection of extraterrestrial intelligence carries no value unless it excludes a plau¬ 
sible model of life in the universe. One might imagine the following model. The Milky Way 
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Galaxy has an age of ~10 Gyr, while only 4.6 Gyr were required to spawn a technological 
species on Earth. Approximately half of the stars in the Galactic disk are older than the 
Sun, and they formed with comparable amounts of heavy elements. Approximately 50% 
of these old disk stars do not have a stellar companion within 100 AU (Duquennoy and 
Mayor 1991), thereby permitting stable planetary orbits. Therefore, ~25% of the stars in 
the Galactic disk have requisite characteristics of adequate age, chemical composition, and 
dynamical quiescence to serve as sites for the development of technological species. The an¬ 
thropocentric nature of these stellar criteria are apparent, and thus serve merely as a useful 
guide for a model. 

Goincidentally, the target stars used here are systematically middle-aged or older and 
void of close companions. They were selected for the Doppler planet search at Keck from 
among stars in the Solar neighborhood. They were chosen to have ages greater than 2 Gyr, as 
judged from the Gall H&K chromospheric emission line, because of the resulting photospheric 
stability of older stars. Moreover the target stars have no known stellar companions within 2 
arcsec (typically ~50 AU) to prevent any companion’s light from entering the spectrometer 
slit (but see Vogt et al. 2002 and Nidever et ah 2002 for a few recently detected, but 
negligibly bright, close companions). Moreover, spectroscopic analyses of the stellar spectra 
show that 90% of the target stars contain amounts of heavy elements that are within a factor 
of 2 of the Solar value ([Fe/H] = -0.3 to -1-0.3; D.A.Fischer, private communication). 

Thus, the majority of the 577 target stars for this SETI search are older than the Sun, 
have roughly Solar chemical composition, and have no perturbing stellar companions within 
50 AU. This stellar sample is clearly biased. Thus our observations test a limited model of 
the Galaxy in which older, single stars harbor technological civilizations. 

The speculative model posits that some fraction of these target stars contain civiliza¬ 
tions that sent probes to, or established colonies on, planets, moons or other platforms either 
in their host planetary system or in other systems. This model of behavior, albeit anthro¬ 
pocentric, stems from the tendencies of Homo Sapiens toward exploration. To communicate 
with such probes or colonies over distances of Astronomical Units to parsecs, some fraction 
of these civilizations might be modeled as using electromagnetic waves, including optical 
lasers. 

We model the solid angle subtended by laser beams as follows. Some fraction, /tech, of the 
total number of target stars, A'"stars=577, harbor technological civilizations. On average each 
such civilization emits some number of laser beams, Nbeams, perhaps in arbitrary directions. 
We consider that each beam subtends an average solid angle Dbeam- 


The total solid angle, Oxotai, subtended by all arbitrarily oriented beams is 
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^Total -^stars/tech-^beams^beam ( 16 ) 

We ignore here details of overlap of the beams, their actual distribution of power and beam 
size, and any purposeful aiming toward or away from us. Any properties of the lasers that 
are constructed with the Earth as a consideration could drastically change the probability 
of interception of such beams. 

In our survey, we have monitored Wtars= 577 stars. The nondetection obtained in the 
survey here suggests that the total solid angle subtended by lasers emitting 50 kW or more 
is constrained as flxotai < 47r sr . 

Thus for arbitrarily directed beams, we constrain the model as 

/techATieamsflbeam ^ 47r/577 (fl^) 

The nondetection reported here suggests lasers emitting detectable power must sat¬ 
isfy the above constraint for the combination of /tech, A^iasers, flbeam- Apparently, advanced 
civilizations communicating by kilowatt optical lasers are not so common as to £11 47r sr. 

As a touchstone, we consider the case of the diffraction-limited laser with an aperture of 
10-m (Keck-to-Keck), giving a solid angle for one beam, hlbeam ~ 10“^® sr. For transmitters 
of such aperture, a laser emitting 50 kW would be detectable. For this case, each star would 
have to emit ~ 10^^ arbitrarily oriented beams on average, each with power P > 50 kW, in 
order that one would likely be oriented, serendipitously, toward us. 

Clearly, the observed nondetection implies that such numerous, narrow beams, however 
implausible, do not exist in sufficient numbers to permit likely detection. These observations 
offer a poor constraint on lasers of arbitrary orientation, as their solid angles are simply too 
small in this model. One might consider a Galactic model consisting of wider laser beams, 
enhancing the chances of interception. In that case, the present nondetection would impose a 
constraint of fewer lasers emanating from the target stars, but would require proportionally 
higher power per laser for detection at G-a above the stellar continuum. 

More meaningfully, the nondetection reported here suggests that none of the 577 target 
stars harbors a civilization that has purposefully directed a laser with P >50 kW toward 
Earth (for the case of a 10-meter diffraction-limited transmitting aperture). This modest 
detection threshold of 50 kW highlights the ease with which an advanced civilization could 
signal us, if desired. But such is apparently not the case. 
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Fig. 1.— Top: A portion of a spectral order from a raw CCD image of the spectrum. Any 
emission line should span the full width of the stellar spectrum, as it is caused by seeing 
at the slit. Bottom: The reduction process converts a two dimensional spectrum into a one 
dimensional spectrum. 
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Fig. 2.— Top: A full chunk of difference spectrum, including spikes, which will be searched 
for laser line candidates. The peak height is 15.80% above the continuum. Bottom: The 
same chunk of difference spectrum with spikes removed and offset for clarity. The standard 
deviation of this spike-free chunk is the reported sigma. 
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Fig. 4.— The observed histogram of the number of segments versus number of hits for the 
star HD29528. The dots indicate the best htting Poisson distribution with the mean number 
of hits per segment equal to 0.00907. 
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Pixels (w.r.t. reference spectrum) 


Fig. 5.— Three hits in the reduced spectra of HD29528. The test spectra (dotted lines) have 
been aligned with the reference spectra (solid lines). Twice the difference spectra are shown 
below. The candidate laser lines are in order 83 and the central peak is in pixel 1436 with 
respect to the reference spectrum. 
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Fig. 6.— An apparent narrow emission line. Top: The raw CCD image of an observation 
of HD29528 taken on JD 2451793.1, centered on the first candidate laser line. Middle: The 
raw CCD image of the reference spectrum. Bottom: The difference of the two raw spectra. 
A bright line spanning the width of the spectral order can easily be seen. 
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Fig. 7.— Top: Order 83 of a B star spectrum. A strong feature can be seen at pixel 1436. 
Bottom: Zooming in on the bad pixels. 
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Fig. 8.— A flaw in the CCD masquerades as an emission line. Top: A chunk of raw 
spectrum from order 83 of the reference spectrum for HD29528 centered on the location of 
the candidate laser lines. Bottom: The same chunk of spectrum with the intensity levels 
adjusted. A dark CCD flaw can be seen under the bright spectral order. 




